set term eps
set output "deuteron.eps"
a=2
V=35
set xlabel "E (MeV)" font "Arial,20"
set ylabel "f(E) (MeV)" font "Arial,20"
unset key
plot [-5:0] 0, sqrt( 1000*(x+V) )/ tan( sqrt(0.025*a*a*(x+V)) ) + sqrt(-1000*x)

